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Abstract 

For a directed polymer in a random medium lying on an infinite cylin- 
der, that is in 1 + 1 dimensions with finite width and periodic boundary 
conditions on the transverse direction, the winding number is simply the 
algebraic number of turns the polymer does around the cylinder. This pa- 
per presents exact expressions of the fiuctuations of this winding number 
due to, first, the thermal noise of the system and, second, the different 
realizations of the disorder in the medium. 

PACS numbers: 05.40.-a, 36.20.-r, 64.60. Cn, 71.55.Jv 



Introduction 

A directed polymer in a random medium is one of the most simple non-trivial dis- 
ordered system and is, as such, of special theoretical importance. Indeed, several 
exact results on directed polymers with strong disorder have been obtained^] 
I2I ISl ^ m 1 so that general approximation schemes developed to tackle more 
complicated disordered systems such as spin glasses could be tested within the 
directed polymer context. The directed polymer is also relevant in the context of 
nonequilibrium phenomena, as it is related, through simple changes of variables, 
to growth models governed by the Kardar-Parisi-Zhang (KPZ) equation[H| E] 
and to non-turbulent flows such as the Astmmetric exclusion process (ASEP) 
model|2lC|. 

The objective of this paper is to study the winding of a directed polymer lying 
on the surface of a cylinder. The algebraic number of turns W the polymer does 
around the cylinder is a random variable which depends on the realization of the 
disorder and which fluctuates because of the thermal noise. The statistics of the 
winding number of a polymer in an homogeneous medium in 2 -I- 1 dimensions 
around a cyHnder goes back to the work of Spitzer|2| and is relevant for the 
physics of vortices in type II superconductors [HI 1101 ITT] . In physical situations, 
the system is, however, usually disordered; the effect of columnar defects has 
been studied analytically ^l], and the winding number around a cylinder of a 
polymer in a random medium with point-like disorder in 2-1-1 dimension has been 
explored numerically (HI E|- When there is an attractive interaction between 
the polymer and the winding center, the polymer can be confined around the 
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cylinder and the system can be regarded as a polymer in the 1 + 1 dimension 
with periodic boundary conditions. In that situation, the present work gives 
exact expressions for the statistics of the winding number. 

The directed polymer on a cylinder is also related to the classical limit of 
strongly correlated fermions in one dimension with disorder (Luttinger Hquids) : 
the X position of the directed polymer corresponds to the phase of the fermions, 
and the phase has, of course, periodic boundary conditions. The winding num- 
ber of the polymer corresponds to the density of fermions. Disorder, while 
periodic in both cases, does not have exactly the same correlations, but the 
models are sufficiently similar to hope for some universalitv [T2l [13] . 

A first result of the present paper states that the thermal fluctuations of 
the winding number are simply equal to what one would obtain for a directed 
polymer in a homogeneous medium; disorder is simply averaged out. A sec- 
ond result concerns the thermal-averaged winding number W. Because of the 
randomness of the medium, this quantity is not zero and the expression of its 
variance {W ) averaged over disorder is obtained. 

The present paper is organized as follows: the first section is a brief recall of 
how the directed polymer in a random medium can be mapped to a quantum 
mechanical problem of interacting bosons using the replica method, and how this 
quantum mechanical problem can be solved with the Bethe Ansatz[Tl ll4l [2ir^fi]. 
In a second section, the winding number is introduced and defined and the two 
main results of this paper are stated in equations l|19ll and l|21|l . Section three 
gives the main lines of the derivation, and, finally, technical points are developed 
in the three appendixes. 



1 Definition, notations and free energy of a di- 
rected polymer 

Let us consider a directed polymer in 1 -I- 1 dimensions where the dimension 
in which the polymer is directed (the "time" dimension) is taken to be very 
large and the transverse dimension (the "space" dimension) has width 1 and 
periodic boundary conditions. As it is directed, the polymer can be described 
by a single- valued function y[t) and the partition function of a directed polymer 
of length t ending at position x is given by 

Z{x,t) = / Vy{s) exp ( - [ds 
Jv(t)=x \ Jo 

where rj^x, t) is the contribution by the random medium to the energy of the sys- 
tem. Disorder in the medium is assumed to be characterized by an uncorrelated 
Gaussian noise of variance 7: 

(77(x,t))=0, {7^{x,t)7^{x\t'))^^5{x-x')5{t-t'), (2) 

where the brackets ( ) represent the average over disorder. 

It is a well known result [21 Hj that this system can be successfully 
mapped to a quantum mechanical problem using the replica method; indeed, if 
we define 

i^{x]t) = il^ixi,. . . ,xn;t) = ^ ^ (3) 
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where 

Z{t) = JdxZ{x,t) (4) 
is the full partition function, then ^ and Q imply 

dip I ^ d'^lp c/ \ , r^x 

i—l * i<j 

with periodic boundary conditions on all the space variables xf. 

'ip{. . . ^ 0, . . .-,1) ^ . . = I, . . .-,1). (6) 

The normalization by (Z)" in (pj is just a simple way to get rid of a low-scale 
divergence introduced by the continuous description (QJ of the system. In other 
words, without this normalization, there would be a trivial extra term in © 
involving the lattice size of an underlying discrete formulation of the problem. 

For an infinitely long polymer, that is, in the large t limit, the amplitude of 
■4>(x;t) is given by the fastest growing mode of JSHEl- In quantum mechanical 
language, we have 

\ogi){x\t) log J dxi .. .dxn^ix;t) 
nm = iim 

t — >oo t t — ^oo t 

^ log{Z{tr)-nlog{Z{t)) (7) 

t^oo t 

where E{n^ 7) is the ground-state energy of the Hamiltonian 

i—l * i<j 

which describes n particles with attractive delta- interactions on a ring of size 1 . 

The same Hamiltonian with a negative value of 7 (that is, with a repulsive 
delta-interaction) has been much studied to determine the spectrum of a gas of 
bosons [T5l IT61IT71IT81 IT9I [2n| . In that context, using the Bethe Ansatz^, it 
was shown that the ground-state energy of ijHJ can be written as 

1 " 

E{n,j)^--J2xl (9) 



2 

a=l 



where the {Aq} are solutions of 

Aq - A/3 -I- 7 



IT 7 7^ . with IimA„=0. (10) 



l<p<n ^ 



Of course, this expression obtained for 7 < remains valid in the directed 
polymer context where 7 > 0. 

In the quantum mechanical problem JHJ, the ground state energy E{n,^) is 
well defined only for integral n; after all, n is the number of particles. However, 
for the directed polymer, (^"), which is related to E{n,'-f) through JJjl, can 
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be defined for arbitrary values of n. The small-rt limit is of special importance 
here: as the directed polymer is a disordered system, the free energy is a random 
variable and (•^") is the generating function of this free energy. Indeed, we have 

!=£<f!>.„0^5f>+!^<i^ + 2!<i^ + o(,,'), (11) 

t t 2 t 6 t 

where ( log (^log^ Z^Jt — (( log^ Z) — ( logZ)^)/t, etc. are the cumulants 
of the free energy per unit-length of the directed polymer. Thus, if we can 
generalize equations ^ IIOII to arbitrary values of n, the expansion of E{n, 7) 
for small n gives, using 0, the distribution of the free energy of the directed 
nolvmer|22|. 

This method was used^ for the directed polymer on a space of infinite width 
in the a;-direction. The Bethe Ansatz equations are then much simpler than 
and one obtains j23l IT], when n is an integer, E(n,j) — 7^(n — n^)/2A. This 
result was used to argue that only the two cumulants ( log Z)/t and ( log'^ Z)^/t 
do not vanish in the large t limit and that, therefore, the fluctuations of logZ 
scale like ^^/^j2ll25ll26l. 

When space has finite width, however, it is easy to see that the free energy 
is an extensive function and that all its cumulants scale like t. In two previous 
papersjHlBI, we solved the Bethe Ansatz equations itTTHl and computed the three 
first terms of the small n expansion of E{n,j). Up to the order n^, the result 
is 

Ein,j)^n{l + l;]^^^ I d\^^^ + 0{n% (12) 
so that, using Q: 



2 24; 4V2 7o tanh^ 

2 V 2 



(■o,Z) ,o,(Z) _/. ^^^^ 



t^oo t V 2 24 



2V2 Jo tanh ^ 

2 V 2 

:j+ll-jL + J^ + 0{"f^) for small 7, ^^^^ 
' 12 360 5040 ^' ' 

^ ^ + 4C(3) + O - for large 7, 
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where C(3) =J2k-^ ~ 1.20206. 



2 Winding number of the directed polymer 

An important topological property of a directed polymer is its winding num- 
ber W, that is the algebraic number of full turns the polymer makes around the 
cylinder on which it lays. One way to define this winding number is to increase 
W by one for each "time" t where the x coordinate of the polymer goes from 1~ 
to 0+ and decrease W by one when x goes from 0+ to 1^. Another way is to 
unroll the x-coordinate and set = / xdt. Of course, the differences between 
those two definitions smear out in the large t limit. 
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As for any quantity in a disordered system at finite temperature, the winding 
number W fluctuates for two distinct reasons. One is the thermal fluctuations: 
for a given realization of the disorder and at finite temperature, the directed 
polymer fluctuates around the path with the lowest energy, and those fluctu- 
ations may change the winding number of the polymer. The other source of 
fluctuations is the quenched disorder on the medium. 

In this work, a horizontal bar is used to denote the thermal average, that 
is the average computed over all the possible directed polymers counted with 
their Boltzmann weights. The cumulants are noted with an extra c subscript: 

W : thermal average of W, 



{W')^ : A;-th thermal cumulant of W, 



with {W^)^ = W^-W , {W^)^ = W3-3WW^ + 2W , etc. These thermal av- 
erages and cumulants are calculated for a given, flxed, realization of the disorder 
and usually depend on that realization. 

The average and cumulants of a quantity Q computed over all the realizations 
of the disorder are written with brackets: 

(Q) : average of Q computed over all realizations of the disorder, 
(Q'')^ : k-th disorder cumulant of Q. 

It is worth noting that, for a given realization of the disorder, the thermal 

average W of the winding number is not zero; the disorder breaks the symmetry 
and may favor one orientation over the other. However, W is an extensive 
quantity and, if we imagine that we cut an extremely long polymer in many very 
long sections, all the sections arc nearly independent and W may be regarded 
as the sum of uncorrelated random variables. Therefore, 

W (W) 

lim — = lim — = 0. (15) 

t— ^oo t t— ^oo t 

This property that W approaches {W) in the large t limit is known as "autoav- 
eraging". Likewise, all the thermal cumulants of W (which are also extensive 
quantities) share the same property: 



Those cumulants, which characterize the thermal fluctuations of a directed poly- 
mer's winding number, depend on the realization of the disorder only when the 
length t of the polymer is finite. 

Other quantities of interest are the disorder cumulants of the thermal-avera- 
ge of the winding number of the polymer. Indeed, the quantity W depends on 
the realization of the disorder, and its fluctuations are characterized by another 
series of cumulants: 

lim ^-A, (17) 

with (W )^ = (W )-{W) , etc. Actually, we might be interested in computing 
many quantities characterizing the winding number, such as 

lim ^ ^ — ^ (18) 
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which represents the fluctuations due to the disorder of the thermal-mean square 
of the winding number, per unit length. 



A first result of the present paper is 

and 



lim 

t — >oo 



t 



1 



lim 

t — ^oo 



t 



for /c 7^ 2. 



(19) 



In other words, thermal fluctuations of the winding numbers are Gaussian and 
independent of the disorder 7. For an infinitely long polymer, the thermal 
fiuctuations of the winding number of the polymer behave as if the directed 
polymer was simply doing a random walk in a disorder- less environment. 

A second result of the present paper is 

t2\ 



lim — lim 

t— ►00 t n— ►O 



-1 



fdE{n,j) 2^, A 1 



- 1 



7 d"f J t 



(20) 



lim 



where E{n, 7) is the ground state energy of the quantum problem computed 
in jni EI and given in l(T2|l . Therefore: 



lim 

t — ^00 



V2 



tanh^ 



+oo ;^4g-AV2 

tanh^^ 



- 1 



7 



8C(3) 



+ 0(7 ^) for large 7, 



(21) 



7 



7 
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0(7^) for small 7. 



(where C(3) = ^ fc-^ ~ 1.20206.) 

The expression of the directed polymer's free energy is written with 
dimensionless variables. If we explicitly put back physical constants and use the 
following expression instead of 

Z{x,t) = jvy{s)e^^iy-P jds ^[^^ + ry(y(s), s) j , (22) 

where (i — (ksT)^^ is the inverse of temperature, k is the rigidity modulus 
of the line and where the spatial dimension x has finite width w and periodic 
boundary conditions, then ifTfljl and l|2Tll become: 



lim 



and 



lim 

t — >oo 



for fc 7^ 2, 



lim 

t — *oo 



(23) 



(3kw' 
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where ^(7) is the scahng function given in II21|I . We obtain the following ex- 
pansions: 

lim — — — — at low temperature, 

« at high temperature. 

360 

3 Derivation of (111) and (EDD 

3.1 Equivalence to a quantum mechanical problem 

To obtain both results ifTfljl and lf2(Hl . we define a new partition function Zz{x, t), 
the purpose of which is to count the winding number of the polymer: 

I — 1 energy of path y{s) 1 -\-z 1 winding number of that patli 1 

Z^{x,t) ^ Vylsje V /V /. (25) 

Jv{t)=x 

The sum is made over all the directed polymers ending in x, and the "energy of 
a path" is the same as in 

Clearly, Zz{t) = J Zz{x, t) dx is related to the winding number W by 



Z,{t)= ZQ{t)e-w. (26) 

If we define the winding number W as an integer that changes by ±1 each 
time the directed polymer wraps around the domain by crossing the a; = 
or X = 1 boundary, then the boundary conditions for Zz is: 

Z,(0,i) = e"Z,(l,t). (27) 

Apart from that, the equations satisfied by Zz{x, t) are the same as the equations 
satisfied by Z{x, t). In particular, if we define 

{ZzAxi,t)...ZzSXn,t)) 
^zi,...,z^[Xl,...,Xn]t) ^ , , n , (28) 

\Za(t)) 

this new wavefunction ip is also a solution of 101; only the boundary conditions 
are changed: instead of ©, the new conditions read: 

^zi , . . (*^1 ; ■ ■ ■ : -^i — 0; ■ ■ ■ 1 "^n 7 ^) — 6 . ..Zn(.-^1 : ■ • ■ : -^i — 1 ; ■ ■ ■ ) -^n : ^) ■ (^9) 

Thus, as in Q, the long "time" t behavior of Zz{t) is given by 
r ^og{Zz^---Zz^)-n\og{Zo) 

hm = -£:(n,7;zi, . . . ,z„), (30) 

t— >oo t 

where E{n, 7; zi, . . . , z„) is the ground state energy of the same Hamiltonian (jsj 
as before, but with the new boundary conditions l|29ll . 

This new ground state energy E contains all the information on the winding 
number W. For instance, from Ij26|l . and by definition of the cumulants, we 
have, for fc > 0, 

W^)c-S^^ogZz{t) . (31) 
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(logZ^) is easily obtained from II3()|I : we set all the {zi} to one single value z, 
make a small n expansion, and retain only the first order. We get: 



(32) 



2 

Getting W is more tricky. We would need (dlogZ^/dz)'^, but that quantity 
can only be obtained from II^OII if the parameters {zi} take at least two different 
values. For example, we have 



/ Zo(l + ziW + 0{zf) 
lim(Z.,Zri)-/ ^ 



\ Zo(l + zT4^ + 0(z2)) / ' (33) 
= 1 + (zi - z){W) - zzi{W^) + 0{zf) + 0(z2), 

and 

limlog(Z,,Zr') - {z,- z){W) - zz,(^{W^)~ (Wy) +0{zf) + 0{z^). (34) 
Therefore, putting all the pieces together. 



hm = hni E{n, 7; zi,z,..., z) 

t^oo t n~*0 dzOZi 



(35) 

z=0 
21=0 



Finally, to obtain the results announced, we need to compute E{n, 7; z, . . . , z) 
and, to the first order in z and zi, E{n, 7; zi, z, . . . , z). 

3.2 Determination of E{n, 7; . . . , z). 

When all the parameters {z^} are equal to one single value z, the problem is 
easy: all the replica play a symmetric role, so that the ground state eigen- 
vector ^z,...,z{xi, ■ ■ ■ ,Xn) of the Hamiltonian (|Hll is a symmetric function of all 
the {xi}. As shown in appendix El the standard Bethe Ansatz derivation gives 
the result. Instead of (Ifll lTTHl . we get 

1 " 

E{n,r,z,...,z)^--Y,^l (36) 



2 

£1=1 

where the {Xa} are solutions of 
If we define 

Aq = Aq + z, (38) 

then the {Ac} are clearly solutions of the standard Bethe Ansatz equations 
Using Ij36ll . we obtain 

E{n,r,z,---,z)^E{n,-f)-^z\ (39) 
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E{n, 7) — E{n, 7; 0, . . . , 0) is the ground state energy Ijl2|l before introduction of 
the {zi}. We have used ^ Aq = 0, which can be easily deduced[3, 4 from lllOII . 

Using Ij32|l , the result on the thermal cumulants of the winding number 
is then immediate. 

This method, based on a Bethe Ansatz, is not the simplest way to ob- 
tain Ijl9ll . Indeed, the result could be obtained using the statistical tilt symme- 
try 27, 28J of the problem; we define the winding number W of a. path y{s) as 
being simply the unrolled coordinate: 



W 



, dy 
as —. 

as 



(40) 



(This new definition is, of course, equivalent to the previous one in the large t 
limit.) The change of variable y{s) = y{s) + zs in the definition ll25|l of gives 
then: 



T>y{s) exp i — jds 



1 / dy \ ^ 



(41) 



Clearly ri{y{s) + zs, s) have the same statistical properties of r](^y{s), s) and one 
gets 



(logZ,(t))-- + (logZo(t)), 



(42) 



from which the result ltT9|l is straightforward. The first derivation with the 
Bethe Ansatz was included here as it demonstrates part of the method used 
to obtain the second result Ij2()|l . which cannot be derived from a statistical tilt 
symmetry argument 



3.3 Determination of E{n, j; zi, z, . . . , z). 

When the parameters {zi} are not identical, the wave function V' is no longer 
a symmetric function of the {2;^} and the problem is much more complicated. 
Therefore, the standard bosonic Bethe Ansatz used in the previous case will not 
work. However, as shown in appendix El using a more general Bethe Ansatz 
that was first introduced to deal with non-bosonic particles (221 EOl E], we get 
the following result: 



1 " 

E{n,r,zi,z,...,z) = 



(43) 



where the {Xa} are solutions of 

Aq - A/3 - 7 



n 

l</3<n 



and where the {Ca} are such that 

r Aq — A/3 



(3=1 



7 



with 



(3=1 



lim Ca — 0. 



lim Xa — 0, 
7— >o 

2^0 



A/3 



(44) 



(45) 
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When zi ~ z, we recover Cq = z, the result of the previous section. 

To determine the fluctuation (W ) of the winding number of the polymer, 
we only need to compute E{n, 7; zi, z, . . . , z) to the second order in the {zi\. 
From H45II . we easily get 

C« = £l±il^ - - zf + 0{{z^Y). (46) 

We deflne, for all a, 

\ J n jw^ 



7 = 7(l-^^b#)- (48) 



and 

7 = 7(^1-- 

771 

Using those new variables into ijiiji . we obtain the famihar Bethe Ansatz equa- 
tions: 

e^°= n ^:'\'^: +oi{z.n m 

l<f3<n Aa - A^; - 7 

SO thatO^l) using the ground state energy E{n,j) given by (jl2ll . 

n n 

^A2 =-2S(n,7) + 0({zj3) and ^ A„ = 0({zj3)_ (gg) 

a— 1 a—1 

From there, using l(T7ll . one can write ^ A„. We flnally get 

1 [zi + {n- l)z]'^ 



E{n,-/;zi,z,...,z) = E{n,-f) 

z n 

"""" ^'\z,-zr + o{{z,y). 



2 n 

1 



t2 



2 ai?(n,7) 
-i;(n,7) 



(51) 



7 ' 97 
Then, finally, from ijSSl, we get the announced result Ij2()ll . 

Conclusion 

Using the replica method with the directed polymer, one obtains a bosonic 
quantum mechanical problem which can be solved by the Bethe Ansatz. By ex- 
tending this method and using a more general Bethe Ansatz that was introduced 
to deal with non-bosonic particles [22l, it has been shown how the different quan- 
tities characterizing the fiuctuations of the directed polymer's winding number 
can be computed using new Bethe Ansatz equations. Building upon a previous 
work[3l B|, those equations were explicitly solved in two cases giving the re- 
sults ifTflll and ISnJHJ ■ The second result is particularly interesting as it simply 
relates through equation lf2?Hl the fluctuations of the thermal-averaged winding 
number and the fluctuations of the free energy of the directed polymer. It would 
be interesting to understand this relation in a more direct way. 
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In principle, the method presented in the present paper should allow to 
compute more cumulants of the winding number and, eventually, its complete 
probability distribution. For that, however, one needs, as a first step, to gen- 
eralize lf5T|l and write the expansion of E{n, 7; zi, z, . . . , z) to higher orders in 
the {zi}. Indeed, one can show that, for example. 



hm = hm E{n,j;zi,z, ...,z) 

t— 00 t n^o azozi 



■ (52) 

z=0 
zi=0 



Obtaining E{n, 7; zi, z, . . . , z) to the fourth order in the {z^} is not, however, 
an easy task, as the trick used in (HTIl would not work at that order. 

As a second step, to compute more complicated cumulants of the wind- 
ing number such as (W ) — 2>(W ) , one needs to generaHze II43H45|I to the 
case where the {z^} take at least four different values. Higher order cumu- 
lants would require, of course, the energy of the system with more different 
values of the {z^}. A matrix approach such as the one developed in the present 
paper could lead to the result. Another possibility might be to try an ap- 
proach similar to the "Nested Bethe Ansatz" method developed by Yang and 
Sutherland (221 EOl 1^ Q] H] to compute the ground-state energy of the system 
described by the Hamiltonian (jsj with different types of particles and symmetry 
relations which depend on the type of particles. Their results are not directly 
applicable to the directed polymer's winding number as all the particles have 
the same symmetry relations but different boundary conditions, but it might be 
worth investigating if the nested Bethe Ansatz could be adapted. 



Acknowledgments 

I would like to thank Pierre Le Doussal, who suggested this work, and Bernard 
Derrida for his interesting discussions. 



A Bethe Ansatz equations when all the zi have 
the same value z 

When all the {z^} are equal to z, all the particles have symmetric roles and 
the standard bosonic Bethe Ansatz leads to the result. To recall the standard 
derivation 1 1 5L [T7] . we look for solutions of the following form: 

V'.,.....(a;i,...,x„) = ^a(a)eSr=iV.)-.w, (Al) 

where the sum is made over the n\ permutations ct of {1, . . . , rt} and where r is 
the permutation defined by 

Xr(l) < Xr(2) < ■ ■ ■ < Xr(n)- (A2) 

The n! amplitudes {a(cr)} and the n pseudo- wave- numbers {A^} are unknown 
variables to be determined. 
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We use this expression of ip in Tiip = Eip, where H is the Hamiltonian JSJ. 
In the regions where all the {xi} are different, it is straightforward to get 



1 " 

E{n,r,z,...,z) = --Y,>^l, (A3) 

a=l 

so that we only need to determine the {Aq.}. At each crossing of two particles, we 
have to ensure the correct discontinuities in the derivatives of to compensate 
for the i5-functions in H. This gives the following conditions, for all a and 
all 1 < A: < n: 

a{a o Tk) = '——a{(j), (A4) 

A(T(fc) - Aa{k+1) + 7 

where Tk is the permutation that swaps k and k + 1 and leaves all the other 
integers unchanged. 

As any permutation a can be written as a product of the elementary per- 
mutations Tk, one can use ljA4ll to write all the {a{a)} up to an arbitrary 
multipHcative factor. However, as the decomposition of a permutation cr as a 
product of Tk is not unique, one must check that the {n — l)n! equations (|A4|I 
are self consistent. The best way to do that is to write down explicitly the 
solution: 

n 'r"'r"'' - m 

l<a</3<n ^ ^ 

It is easily checked that this is indeed the solution of all the equations IIA4|1 . 

So, for any set of values {Xa}, the wave function ijAljl where the {a((T)} are 
given by l|A5ll is an eigenvector of the Hamiltonian jHJ. The values of the {Aq} 
can then be obtained from the boundary conditions l|29ll . One gets: 

a(cr) = e^+^"(i)a(croC), (A6) 

where C is the circular permutation C(l) = 2, C(2) = 3,. . . , C{n — 1) — n, 
C{n) — 1. Using ljA5ll . we easily get the new Bethe Ansatz equations. For all a, 

e^°+^= n + (A7) 

We are only interested in the ground-state solution. By continuity of this 
ground state, we get the last condition 

lim Aq = 0. (A8) 

B Bethe Ansatz equations when all the {zi\ ex- 
cept zi have the same value z 

When the parameters {zi} do not take the same value z, the computation of 
the energy E is more complicated; indeed the wave function -i/i is no longer a 
symmetric function of the {xi\ and there is no way that the standard Bethe 
Ansatz ijAlll might lead to the result. 
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However, in order to study the Hamiltonian Q for fermionic particles or, 
more generally, for particles with arbitrary symmetries and anti-symmetries, 
a more general Ansatz than llAl|l has been proposed[22| I30L 1^: in llAl|l 
and ljA2p . the permutation r is only introduced as a convenient way to get 
the coordinates {xi} of the n particles sorted from the leftmost particle to the 
rightmost in the expression of the wave function. An easy way to break the 
symmetry of is to make the parameters {a{o')} explicitly dependent on the 
permutation r: 



5 ) 



(Bl) 



where the permutation t is, as before, defined by ljA2|l . 

As shown below, the solution to our problem with the unusual boundary 
conditions Ij29|l can also be written using IIBlll . We first begin with the most 
general case where all the {zi} are different and, at some point, specialize to the 
simpler case where all the {zi} except zi are identical. 

B.l General Bethe Ansatz equations for arbitrary {zi}. 

Using the Ansatz (|Blll in H^j = Eip where H. the Hamiltonian (jHJ, we have, as 
usual: 



E{n,'-f;zi, 



(B2) 



The new equations for the parameters {a(T, cr)} are more complicated than ljA4ll : 

(ACT(fe) - Xa(k+i))a{T o Tk, cr) - -yair, a) 



a(T,ao Tk 



(B3) 



T(fe) - K{k+i) + 7 
for any permutations t and a and for any integer 1 < k < n. 

A convenient way to write the (n!)^ parameters {a(r, cr)} is using nl vectors 
indexed by a, each vector having n\ components: 



a(CT) 



a(ri,cr) 
a(T2,cr) 

a{T„\,<j) 



(B4) 



where ti, . . . , t„! are the nl permutations of {1, . . . , n} sorted in an arbitrary 
way chosen once for all. (The order must, of course, be the same for all values 
of a.) We now introduce the matrices Mk defined by 



a{Ti o Tfe, (t) 
a{T2 o Tfe, (t) 

a(T„! oTk,cr) 



= Mka{a). 



(B5) 



Those matrices just shuffle the components of the vector d{a)\ there is thus 
exactly one "1" per raw and per column and all the other components are "0". 
In a concise way, we can write Mk as 



(Mk),^, 



51'. 



(B6) 



13 



The matrices are a representation of the permutations Tk- As such, they 
have the same standard commutation properties as the permutations: 



Ml = /, MkMk+iMk = Mk+iMkMk+u 



(B7) 



MkMk'=Mk'Mk if|fc-A:'|>l. 
{/ being the identity matrix.) Equation ljB3|l is then simply written as 

a(aor,)=if^'^'"^'+''aM, (B8) 
where Yj^'-' is the Yang-Baxter operator defined |29] by 

It is clear from IIB8II that any vector a{a) can be obtained from the knowledge of 
one of them. However, as in the symmetric case, one has to check that the result 
does not depend on the way the permutations are decomposed as a product of 
the elementary permutations Tk . There are no explicit formula 32; such as (|A5|I 
of a(cr), but one can check that the {n — l)n\ relations ljB8|l are indeed self 
compatible. This is implied by the following "Yang Baxter" relations j29j 

Yk'^S'"' =Y^/''y^^' if|A:-fc'|>l, (BIO) 

^k ^k+l^k ~ ^k+l^k -'fc+1' 

which can be easily checked using IIB7|l and ijBQp . With the first of those three 
relations, using ljB8|l twice to compute a(cr o Tk o Tk) gives correctly a{a). The 
second relation impHes a{aoTkoTk') = a{aoTk' oTk) if \k — k'\ > 1 and, finally, 
the third relation gives a{aoTkoTk+ioTk) — a{aoTk+ioTkoTk+i). It is a well 
known property of the symmetric group that those three necessary conditions 
are actually sufficient to ensure that the relations (|B8|I are self consistent. 

One still needs to write the boundary conditions Il29p with the parame- 
ters {a{T,a)}. From IIBlll . one gets 

a(T, a) — exp(zT.(i) + Acr(i))a(T oC,a o C). (Bll) 

where C is, as in IIA6|l . the circular permutation. 

As C = T10T20 ■ ■ -oTn-i, the matrix that shuffles the lines of the vectors a (a) 
according to the permutation C is simply the product of the matrices Mk. Thus, 
we have 

a(a) = exp(A<,(i))ZMiM2 . . . M„_ia((T o C), {B12) 
where Z is the diagonal matrix defined by 

(Z),,, =5^exp(z,.(i)). (B13) 

Moreover, using several times IIB8|l . we get, from the definition of C, 

d{a o C) = y-W-Wy-W -("-!) _ y-(i).-(2)5(^). (B14) 
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Putting together IIB12|l and ljB14ll . we see that a,{a) must be, for each tr, the 
eigenvector of some operator. 

exp(A,(i))ZMi . . . M„_iy,r«'"(") . . . Y^'^^^^^^^^ia) = a{a). (B15) 

There exists a non-zero a{a) such as ljB15|l holds only for certain values of 
the {Aq}. However, as the n\ vectors {a((T)} are not independent variables, 
we need to check that the n\ relations ljB15|l are compatible: they must hold 
simultaneously for the same values of the {Aq,}. As explained in appendix El 
this is the case. To obtain the values of the {Aq} we are looking for, we use 
when a is the identical permutation. A non-zero a(cr) exists if and only if 

det (/ - exp(Ai)ZMi . . . Mn-iY^-^\Y^,-^^^ . . . y^^'') = 0, (B16) 
or, using the properties l|Bir)|l of the operators F, 

det (y^^Y^^^ . . . y„"i\ - exp(Ai)ZAfi . . . M„_i) = 0, (BIT) 
or, using the definition l)B9|l of the operators Y, 

pj^^ (Ai -A.+i)Mq+7/ _ ^ A ^ 0. (B18) 

\a—l a—1 / 

That last equation relates exp(Ai) to the {Aq}. There are n—1 other equations 
giving all the exp(A/c) which we can obtain either by using ljB15p with differ- 
ent permutations a either, as they play symmetric roles, by shuffling the {Aq} 
in ||bT81 . 

Finally, the wave function ijBip introduced is indeed an eigenvector of the 
Hamiltonian lO with the boundary conditions Ij29|l . provided that the {Aq} are 
such that ijBlSIl and the n — 1 other relations obtained by symmetry hold. 

Note that the new Bethe Ansatz equation ijBlSII can be regarded as a poly- 
nomial of degree n! in exp(Ai), so that we have not one value of exp(Ai) as a 
function of the {Aq}, but n\. This could be expected, as the method we have 
used is known to generate, when Zi — 0, not only the bosonic solution, but 
all the eigenvalues of ((HJ for arbitrary symmetry relations between the parti- 
cles. So, if we write from IIB18|l the n\ possible expressions of exp(Ai) and make 
the {zi} go to zero, we will recover the usual bosonic Bethe Ansatz solution ltTn)l . 
but also the fermionic solution exp(AQ) = 1 and all the intermediate cases. In 
our problem, we are looking for the ground state energy of discernible particles 
and it is known, in this situation, that the ground state is given by the bosonic 
solution. 

To sum up, what remains to be done is to single out from ijBlSII the ex- 
pression of exp(Ai) which goes to the standard bosonic equations lfT?Hl when 
the {zi} vanish, to write by symmetry the n—1 remaining equations giving all 
the {exp(Aa)} as functions of the {Aq}, to solve those non-algebraic equations 
in order to write the ground state energy E — — (1/2) -^q; 3Jid, finally, to 
take the limit n ^ and various derivatives with respect to the Zi to obtain the 
different quantities characterizing the winding number of the polymer. 

As this seems to be a difficult task in the general case, we will go through 
this program only when all the {zi} have the same value z except for zi. 
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B.2 Simplification when all the {zi} are equal except for Zi. 

When all the {zi} are set to zero, the matrix Z is the identity matrix and one of 
the solutions of IIB18|I must be the standard bosonic Bethe Ansatz equation (tlllll . 
One way to see it is to notice that, when ifTTHl holds, the vector which cancels 
the matrix in ljB18|l is simply (1, . . . , 1). Another way to see it is to notice 
that to derive ijBlSII . we never actually used the matrix representation of Mk, 
yj!'"', etc., but only the commutation properties of those matrices. If we were 
to choose other representations of those matrices having the same commutation 
properties, relation ijBlSII would still be valid with those representations. For 
example, when all the {zi} are zero, the bosonic solution is obtained from IIB18|l 
by choosing the trivial representation Mfc = 1. The fermionic solution exp(AQ) — 
1 is obtained by choosing Mk = — 1, etc. 

In the situation where Zi = z for i > 2, with only zi different from z, we 
can make a similar simplification. Indeed, in that case, the particles X2, ■ ■ ■ , Xn 
play symmetric roles. Thus, the ground state solution must be symmetric in 
those variables. Back to the wave function ljBl|l . this means that the parame- 
ters a{Ti,a) and a{T2,<T) must be equal if r]~^(l) = r^^(l). In other words, the 
parameters air, a) do not depend on the whole shuffling r of the particles {xi}, 
but only on the position of xi relatively to the other; for each a, there is one 
value of a(T, a) corresponding to the first particle being the leftmost, another 
value when the first particle is the second leftmost, etc. 

This suggests that equation IIB18|l can be written in that situation with a 
representation of the Mk as matrices of size n x n instead of n\ x nl. Indeed, 
we write the new vector d{a) as 



a(T, cr) such that t(1) = 1 {xi is the leftmost particle) 
a(T, cr) such that t(2) = 1 {xi is the second leftmost particle) 

d(a] 

a{T, cr) such that T{n) = 1 (xi is the rightmost particle) 

(B19) 

Then, the matrix Mk that switches the fc-th and k + 1-th particles is given by 

pl+i ifi-fc, 
{Mk)^^J = <S'; + (B20) 

I (5f otherwise, 



that is: 



Ml 



/O 
1 



10 


10 
1 



M. 



/I 

10 

10 

1 

V 



(B21) 



etc. With all the {zi} but zi equal to z, we can also write the new matrix Z in 
this representation. It is a diagonal matrix, on the first line we have exp(2;i) as 
xi is then the leftmost particle, and, on the other lines, we do not know which 
particle is the leftmost, but it is of no importance as we know it is not xi and 
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as all the other particles have the same parameter z. Thus, we have: 



Z = 



/e^i 




V 



(B22) 



0/ course, the new matrices Mk have the correct commutation relations ljB7|l 
and the final result l|B18|l is still valid with the new matrices Mk and Z. 



B.3 Explicit expression of the determinant 

When Zi — z for i > 2, using the new matrices Mk and Z, we can compute the 
determinant in ljB18|l by induction. First, we normalize exp(Ai) by the standard 
Bethe Ansatz expression and write 



=Ai+Ci _ TT ^1 - +7 



(B23) 



where Ci is the quantity we are trying to determine. 

Using this new variable, we write the determinant in (IB18|I . up to a multi- 
plicative prefactor, as 



dn = det (An - exp(-Ci)-B„), 



where 



, _ TT (Ai - Xa+i)Ma +7 
""11 A1-A.+ 1+7 



(B24) 
(B25) 



and, as can easily be seen. 





/o 


. 


. 




n-l 




. 


. 





B^ = Z\l M,, = 





. 


. 







[0 


. 




0/ 



(B26) 



As appears only once in the matrix, the determinant d„ can be written as 

dn = an-e'^Pn, (B27) 

where a„ and /?„ do not depend on zi. 

Going from An to An+i is easy enough; singling out the last term in the 
product IIB25|l . one has 



An+l 




\ 



(0) 



(0) 



1 ~ Mri+l A*n+1 / 



(B28) 
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with 

^^n+l = T J ^ ■ (B29) 

Ai - A„+i + 7 

Doing this last multiplication, we see that the first n — 1 columns of A^+i are the 
first columns of A„ padded with one final zero, and that the n-th and n + 1-th 
columns of An+i are the n-th column of An with different multipHcative factors 
(respectively fin+i and 1 — ^J-n+i) and different paddings (respectively 1 — 
and /in+i.) Thus, if we develop dn+i over the last line, there are only two terms 
which look very much like d„. The only differences are that the last column 
is multiplied by some factor and that the term e^^ is either missing or not 
multipHed by the numerical factor that affects its column. Finally, using ljB27p . 
one can get 

dn+i = ^„+i(^„+ia„) - (1 - ^„+i - e^"^i)((l - fj.n+i)an - /3ne^'), (B30) 

or 

Un+l = (2^„+i - 1 + (1 - ^J,n+l)e^~^^)an, Pn+l = "(l - ^J.n+l - e^~'^i)/?„ 

(B31) 

with ai = 1 and /3i = exp(— Ci). It is now easy to compute a„ and /?„: 

_W (Ai-A,)(e-C^-l)-H7 

"""U A1-A.-I-7 

'-^ (B32) 

fi Ai - A. + T 

and, finally, the condition that the determinant ijBlSII is zero gives the following 
result: 



71 r 

i=l 



Yl i + ^i^(e-0_i) ^g.iJI e^-«^ + ^^i^ (e^-^^ - 1) 



Ai — \i 



(B33) 

As all the {Aq} play symmetric roles, this is exactly the result announced lUfijl . 

C Proof equation ( IB 151 ) can be satisfied simulta- 
neously for all permutations a. 

To prove that equations ljB15|l are indeed compatible, we start by assuming that 
the {a(cr)} are such that ljB8|i holds for any a and k. As a consequence, 
is true for any permutation and ijBlfip is equivalent to ljB12ll . 

Furthermore, we assume that ljB15ll (or (IB12|l ) is true for a given permuta- 
tion a. To show that it is also true for any other permutation, it is sufficient, 
by induction, to prove that (|B15|I (or (|B12|l l holds for a oTk with 1 < k < n. 

It is necessary to distinguish the two cases k 1 and fc = 1: 

• When 1 < k < n 

When k ^ 1, we have the following properties: 



TkoC = CoTk-i, y;:'z = zyI'' 



(CI) 

Y^''MiM2 . . . M„-i = M1M2 . . . M„_ii;!f 
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The first relation is a basic property of permutations, the second relation comes 
from the fact that Y^'-' does not change the value of t(1), and the third one, 
considering the definition ljB9|l of Yfc , is a rewriting of the first relation in the 
matrix representation. 

We can now show that IIB12|I and, therefore, l|B15|l holds for a oTk'. 

= y-W'-C^+i) exp(A,(i))ZMiM2 . . . M„_ia(a o C), 

= exp(A,(i))ZMiAf2 . . . M„_iy,"J'[)'"("+'^a(a o C), (C2) 

= exp{K(i))Z M1M2 ■ ■ ■ M„_ia(CT o C o Tk-i) 

= exp(A<,(i))ZMiAf2 . . . Mn-id{a o Tfc o C). 

As (t(1) = (7oTk{l) for fc 7^ 1, this is indeed ljB12ll applied to the permutation cro 
Tfc. 

• When k = 1 

Equation ijBlfill express that a{a) is an eigenvector of 

A = ZAh ■ ■ ■ M„_ii;riV'"("V„"iy . . . y-(i)-(3)y-(i)'-(2)^ (C3) 
and we want to prove that d{a oTi) is also an eigenvector of 

ZMi M ^^'^(^'''^("^y'^^^^''^^"^^) -yA'T(2),f(3)yCT(2),iT(i) ^^^^ 

As a(cr o Ti) = y]°^*'^-'''^'^^'a(f7), this is equivalent to prove that d{fj) is an eigen- 
vector of 

B = ri"('^^"(^)zMi . . . M„_iy„"if '"("V„"i'2)'"("-^) . . . 12"''^'"^^^ (C5) 

(The relation Y^'-'Y^''^ = I has been used twice.) 

To conclude, we presently show that AB = BA, which implies that A et -B 
have the same eigenvectors. 

First, we define another diagonal matrix Z2 by 

(Z2),^^=<5^xp(z.^(2)). (C6) 

(Compare with ijB 1311 .1 Clearly, we have 

MiZ = Z2M1 and M1Z2 = ZMi. (C7) 

Moreover, as a consequence, the product ZZ2 = Z2Z commutes with Mi and Yi. 

When computing AB, two matrices Yi cancel out. The matrix Z commutes 
with all the the Yfe and all the Mk except Yi and Mi, so that we can "move" 
the second Z to the left and obtain 

AB = ZiH M,)Y"iV'"*"^ . . . Kr''''"*^^^(n ^^^)5Cii^'"^"^ • • ■ ^2"^"^"^'^ : 

= ZZ2il[ A^OC-V^"^"^ • • ■ Kff^^^'^^) ([[ Af,)y„"if . . . Y^"^'^'"^^^ . 

(C8) 
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When computing BA, the second matrix Z can also travel to the left; we get: 



X . . . J-j 



(C9) 



Thus, the products AB and Svl share the same prefactor ZZ2, so that if 
AB = BA is true when Zi = (the case studied by Yangj^Hl), then AB = BA 
is true for arbitrary values of the {zi}. As it is a well known fact that the 
operators commute in Yang's case, we could stop the proof here. However, for 
completeness, let us properly finish it. 

We continue the simplification of AB; using we can have the whole 

first group of matrices in AB go to the right through the second prod- 
uct Ml . ..Mn-i- We get 

AB = ZZ2(n M,)'r,r-2 '"'"^ ■ ■ • r;(l)^^(3)y.^2),.(n) _ ^^^(2),.(3)_ ^^^^^ 

We do the same for the product BA: 

(Cll) 

Using 

r/^^(Mi . . . M„^i)2 = (Afi . . . Af„_i)2y,::^;, (ci2) 

which can be deduced from the properties ljB7|l of the matrices Mk, we get 

(C13) 

AB and BA have the same prefactorZZ2(n -^^i)^; we need to show that the 
two products of matrices Yfc are equal. We proceed by induction: It is clear 
for n — 1 (or n — 2) and we assume it is true for n — 1. In both products, we 
"move" the matrices l^-i to the left. We get: 



^ ^n-2 ■ ■ ■ ^2 



2^(T(l).<T(n)^o-(2).<T(n)^a(l),fT(n-l) ^a{l),a{3) 

(C14) 

^cr(2),o-(n-l) ^(T(2),cr(3) 
r 

and, using ijBl 

BA = zZ2(J^Mi)'^Y''^'f''^''^''Y'^^^^''^^'^''Y'^^^-^''^^"^Y'^^'^^''^^"'^'' y|^(2),o-(3) 

^ff(l),<T(n-l) ^<t(1),£7(2) 
^ ^n-2 . . . -Tj , 

= ^^2(]^Mi)^y'^*y''^'-"V'^*^i'''^'-"V'^'-2^''^^^^y°"''3^''^^""^-' y^'^(2),o-(3) 

^£7(l),<T(n-l) ^a(l),a(2) 
X i-„_2 . . . -Tj 

(C15) 

Leaving aside the common prefix ZZ2{Yl Mi)"^, the products AB and BA start 
with the same two Y matrices, and what remains are the products of Y matrices 
in the expressions of AB and BA at order n — 1. This, by induction, proves 
that AB = BA. 
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Finally, putting everything together, we have shown that the n\ proper- 
ties ljB12|l obtained from the boundary conditions are self compatible. 
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